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ABSTRACT 

Recent work has linked the quiescent luminosities and hard/soft spec- 
tral state switches of neutron stars (NSs) to their spinning magnetic 
fields. It is shown here that the quiescent luminosities and spectral 
state switches of galactic black hole candidates (BHCs) could be pro- 
duced in the same way for spin rates below 100 Hz and magnetic fields 
above 10^^ G. It is also shown that the ultrasoft peaks and large fiick- 
ering amplitudes of the BHCs would be expected from the surfaces of 
massive NSs. None of the few spectral characteristics that distinguish 
BHCs from low mass NSs have been explained in terms of event hori- 
zons. Serious consideration of the possibility that they might simply be 
massive NSs opens an avenue for proof of event horizons by negation, 
but requires the use of a space-time metric that has no event horizon. 
The Yilmaz exponential metric used here is shown to have an inner- 
most marginally stable orbit with radius, binding energy and Keplerian 
frequency that are within a few percent of the same quantities for the 
Schwarzschild metric. A maximum NS mass of ~ lOM© is found for 
the Yilmaz metric. The two metrics essentially differ only by the pres- 
ence/absence of a surface for the BHCs, thus enabling proof or disproof 
of the existence of event horizons. 

Subject headings: Accretion, Black Hole Physics, Stars: neutron. Stars: 
novae. X-rays: stars 
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1. Introduction 

According to General Relativity, an ob- 
ject of nuclear density and more than 2.8 
M0 would be a black hole (Kalogera & Baym 
1996, Friedman & Ipser 1987). Several black 
hole candidates (BHCs) have been found in 
x-ray binary systems. At this time there 
are ten galactic x-ray sources known to ex- 
ceed the neutron star (NS) mass limit. Only 
two of these, Cyg X-1 and LMC X-3, are 
not in low mass binary systems. Perhaps 
twenty more candidates have been identified 
via spectral similarities (e.g. see Barret, Mc- 
Clintock & Grindlay 1996). Active galactic 
nuclei (AGNs) are also generally believed to 
exceed a Schwarzschild mass limit, however, 
they will be considered here only in passing 
in order to restrict the scope of a lengthy ar- 
ticle. In addition, their radii are less well 
constrained than those of the x-ray novae. 
Nevertheless, it should be noted that there 
are strong spectral similarities between NSs, 
BHCs and AGNs. Campana et al. (1998a), 
Becker & Triimper (1998) and van der Klis 
(1994) have recently reviewed NS properties, 
while Poutanen (1998) has reviewed BHCs 
and AGNs. 

X-ray novae display a rich mix of spec- 
tral and timing characteristics. Many of them 
are keyed to the luminosity level. Luminos- 
ity generally derives from accretion of mass 
from a companion star via an accretion disk. 
Flares are caused by instabilities of the flow 
either from the companion or through the 
disk. Maximum luminosities during major 
flares are often ~ 10^ — 10^ times quiescent lu- 
minosities. Maximum luminosities are often 
near the Eddington limit (~ 2x10^^ erg/s for 
a canonical 1.4 Mq NS). Complex changes of 
luminosities in soft x-rays (0.1-4.0 keV), hard 
x-rays (4-20 keV) and hard tails, (20-200 keV) 



occur during flares. Hard tail luminosity with 

an inverse power law dependence on photon 
energy is prominent during the initial rise of 
luminosity and also during the decay. In most 
cases this power law component weakens sub- 
stantially at high overall luminosities. "Ul- 
trasoft" radiation (White & Marshall 1984) 
with a general brehmsstrahlung shape and a 
peak near 1-3 kev is often, but not always, 
seen at high luminosities. Observations dur- 
ing "dips" where some sources, seen at large 
inclination through absorbing materials, pro- 
vide some spatial resolution, show that the 
soft peaks arise near the center of the the 
disk. BHCs simultaneously show hard tails 
and softer, more luminous peaks more com- 
monly than NSs. Spectral state switches oc- 
cur between soft "high" luminosity states and 
the harder "low" states that are character- 
ized by power law spectral features. Hard 
low states typically have luminosities below 
]^q36 Qj.g^g ]sjgg g^j-^^j below 10^'' erg/s for 

BHCs. Campana et al. (1998) and Zhang, 
Yu & Zhang (1998) have attributed the spec- 
tral state switch to "propeller effects" (lUar- 
ianov & Sunyaev 1975) of magnetic fields for 
NSs. Black hole models usually attribute the 
switch to accretion disk instabilities. In either 
case, it is generally accepted that the hard 
spectrum of the low state originates in the ac- 
cretion disk and that it persists, usually be- 
coming harder, as quiescence is approached. 
There is no consensus about how it is pro- 
duced. The strong similarities of quiescent 
NSs and BHCs (Tanaka & Shibazaki 1996) 
suggest that a common mechanism produces 
the quiescent luminosity. 

Low frequency, ^ 6 Hz, quasi-periodic os- 
cillations (QPOs) often appear in suitably 
high states (Fortner, Lamb & Miller 1989, 
Miyamoto, Kimura & Kitamoto 1991, Mak- 



2 



ishima et al. 1986) of both NSs and BHCs. 
Large amplitude flickering 10~"^ - 10 s lu- 
minosity variation) occurs when power law 
components are dominant, particularly at in- 
termediate luminosity (Balucinska-Church et 
al. 1997, Tennant, Fabian & Shafer 1986, 
Stella et al. 1985). The variations of hard and 
soft x-rays are often correlated (usually with 
delayed hard photons) on short time scales 
(< Is), but are uncorrelated (Pan et al. 1995) 
or anti-correlated on scales of hours or more. 
All of these features, some of which were once 
thought to be black hole signatures, have been 
observed in both BHCs and NSs (e.g., see 
Singh et al. 1994, Barret, et al. 1992, Chu- 
razov, et al. 1995, Stella et al. 1985, Tanaka 
1989, van der Khs 1994). Cir X-1, which was 
dropped from the BHC list after displaying 
surface thermonuclear bursts, shows most of 
these characteristics. Its wide array of spec- 
tral characteristics is likely due to a rather 
wide variation of accretion rates associated 
with an eccentric orbit. 

Some NSs show some features that are not 
obviously shared by BHCs. Radio or x-ray 
pulsations, and surface thermonuclear bursts 
are NS signatures that have only been found 
in erstwhile BHCs such as Cir X-1. NSs usu- 
ally show bursts or pulses but not both. The 
accepted explanation for this dichotomy is a 
magnetic field, a NS signature, that can con- 
centrate accreting matter in magnetic pole 
regions where fusion proceeds at rates high 
enough to suppress bursts (Taam & Picklum 
1978). Bursts are generally not observed for 
pulsars with ~ 10^^ G magnetic fields, but 
have been observed for millisecond pulsars, 
which have weaker fields. Interaction of the 
accretion flow with the NS magnetic field is 
believed to account for the characteristics of 
color-color diagrams and hardness-intensity 



diagrams of "Atoll" and "Z" type NSs. These 
sources have magnetic fields of ~ 10^ G and 
10^ G, respectively (White and Zhang 1997). 
The peak accretion rates of Zs are near Ed- 
dington limits, while those of Atolls are much 
lower. Van der Klis (1994), who developed 
much of the Atoll- Z phenomenology has pro- 
vided a review and a discussion of the sim- 
ilarities between the various states of BHCs 
and NSs; particularly the strikingly similar 
time-resolved power density spectra of flick- 
ering Atolls and BHCs. 

The presence of a hard spectral tail at over- 
all luminosity levels above 10^^ erg/s is con- 
sidered to be a reliable signature of a BHC 
(Barret, McClintock & Grindlay 1996). Why 
an accreting black hole should have this prop- 
erty is not understood (Kusunose, Minishige 
& Yamada 1996, S.N. Zhang, et al. 1997). 
With this exception, the lack of spectral or 
timing signatures of black holes has left mass 
determination as the only way to certify a 
BHC (McClintock 1998). Until signatures of 
event horizons are found, it remains feasible 
to suppose that the BHCs might simply be 
massive neutron stars. The point of this work 
is to explore this possibility. 

BHCs as compact, massive NSs could pos- 
sess magnetic fields, and should have large 
surface binding energies and substantial red- 
shifts of surface radiations. Magnetic fields 
would permit a common mechanism for the 
spectral state switches of BHCs and NSs. 
Larger surface binding energies and redshifts 
for the more massive BHCs would make their 
surface emissions both brighter and softer. It 
would seem necessary to reject these possibil- 
ities before acceptance of the reality of event 
horizons for the BHCs. In order to pursue 
these possibilities in a quantitative way, a 
space-time metric that has no event horizon 
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is needed. The Yilmaz metric (YM) has been 
adopted here for this purpose. The massive 
NS hypothesis discussed here is only weakly 
dependent on metric properties in general, 
but depends crucially on the distinctions be- 
tween surfaces and magnetic fields vs event 
horizons. It is shown here that YM and SM 
accretion disks have such similar properties 
that they essentially differ only by the pres- 
ence/absence of a surface. 

Although most of the mathematical details 
have been placed in an appendix, a compar- 
ison of the two metrics is given in Section 2. 
Results for the two metrics for a simple model 
of a star of constant proper density (Clapp 
1973) are also given. An object of this sort 
with a density in its outer layers correspond- 
ing to nuclear saturation density provides a 
rough approximation of a neutron star with 
a fairly stiff equation of state. It provides 
very conservative estimates of the radiant en- 
ergy and redshifts to be expected for accretion 
reaching a massive NS surface. Magnetic field 
phenomena that affect the accretion flow and 
spectral characteristic are outlined in Section 
3. Section 4 discusses magnetic and strong- 
field gravitational effects in NSs and BHCs. 
Section 5 considers the origins of various spec- 
tral characteristics and the implications for 
massive NS models. 

2. The Yilmaz Metric 

While not widely known, the Yilmaz the- 
ory (Yilmaz 1958, 1971, 1975, 1992, 1994, 
1995) has some interesting positive features 
that may help to resolve the hiatus between 
General Relativity and Quantum Mechanics. 
At this time, it is quite clear that one or both 
of these landmark theories of the twentieth 
century must be modified before they can be 
fully compatible. The Yilmaz theory modifies 



General Relativity primarily by the explicit 
inclusion of the stress-energy of the gravita- 
tional field as a source of space-time curva- 
ture. With a true field stress-energy the Yil- 
maz theory possesses a field Lagrangian and 
can be quantized (Yilmaz 1995). 

Although the Yilmaz theory has been crit- 
icized, (Will 1981) the criticism appears to 
be based on an incorrect assumed form of 
the metric (Yilmaz 1981). For the present 
purposes, the essential features of any alter- 
native metric would be lack of event hori- 
zons and readily calculable NS characteristics. 
Nevertheless, the Yilmaz theory passes the 
four classic weak-field tests. It permits local 
energy-momentum conservation and has no 
adjustable parameters, no singularities and 
no event horizons (Alley 1995, Yilmaz 1994). 
Gravitationally compact objects can exist in 
the YM but they are not black holes. Radially 
directed photons can always escape. 

The metric of space-time in the vicinity of 
neutron stars is dominated by the gravity of 
the star which is considered here to be a static 
object. The static limit interval in the Yilmaz 
exponential metric is: 

ds^ = g{rydf-{dr''+r^de^+r^sm^ed(l)'')/g{r) 

(1) 

where: 

g{r) = exp(-2ii(r)) (2) 

and u(r) is the gravitational potential. In the 
SM: 

ds^ = g{rydf - dr^gir) - r^dO^ - sin^ ^^0^ 

(3) 

and: 

g{r) = 1 - 2u{r) (4) 

The gravitational potential at distance r from 
mass M, u{r) — GM/c^r reaches ~ 2xl0~^ at 
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the photosphere of the sun, ~ 3xl0~^ in the 
solar system relative to the Great Attractor, 
~ 10~^ at the smallest AGN radii of which 
we can be certain, ~ 0.17 at the innermost 
marginally stable orbit of a NS or BHC, ~ 
0.25 at the surface of a maximum mass NS in 
the SM and 0.5 at the event horizon of a black 
hole. The five orders of magnitude jump from 
the solar system testing ground to the neutron 
star regime is good reason for caution about 
the acceptance of any strong-field theory of 
gravitation. 

The exponential metric can be inferred 
from special relativity and the principle of 
equivalence applied to frames co-moving with 
an accelerated particle (Einstein 1907, Rindler 
1969, Yilmaz 1975). This is a strong indica- 
tion that the exponential metric is the proper 
static limit for gravitation. That this metric 
does not satisfy the field equations of Gen- 
eral Relativity was one of the motivations for 
the earliest version of the Yilmaz theory (Yil- 
maz 1958). It is noteworthy that the the ex- 
ponential metric is merely rcscaled by addi- 
tion of an arbitrary constant to the potential 
whereas the SM depends on an absolute po- 
tential. Without this dependence, event hori- 
zons could be transformed away. With the 
conventional choice of zero potential at infin- 
ity, however, the two metrics are the same to 
first order in u(r). This is sufficient to insure 
that they give the same results in the four 
classic weak-field tests of General Relativity. 
Gravitational redshift of radiation observed 
distantly is given by (Rindler 1969): 



z — exp{u{r)) — 1 
and in the SM: 

1 



- Mr) 



(5) 



(6) 



(i) Accretion Mechanics 

In either metric there is an innermost marginally 
stable orbit that can be reached by an ac- 
creting particle. The orbital radius is (see 
appendix) Vms = 4/(3 - \f^)GMl(? in the 
YM. Thus u{r^,) = 0.191. In the SM r^^ = 
6GM /c^ and u{rjn,^) = 1/6. An accreting par- 
ticle starts from a large radial distance, with 
essentially zero momentum and u{r) ~ 0. Its 
energy is E — rrioc^. At r^s in the YM, 
the energy is £■ = 0.945moc^ (see appendix). 
The binding energy at Tms is the difference, 
0.055moc^. Whether this is radiated or ad- 
vected depends on the opacity of the accre- 
tion disk and the viscous dissipation mecha- 
nisms. For any circular orbit in the disk, the 
particle energy in the YM can be shown to 
he E = moc2exp(-M){(l - m)/(1 - 2u)}^/\ 
Thus, in the YM the maximum fraction of 
accretion mass energy which can be observed 
as radiation from the disk is: 



I 1-u 
l-2u 



(7) 



limited to u < 0.191. For the SM the maxi- 
mum fraction of accretion mass energy which 
can be radiated from the disk, limited to u < 
1/6 and 0.057moc2, is: 



fd 



1 - 



2u 



3u 



(8) 



If the innermost marginally stable orbit is the 
inner disk boundary, disk mechanics and disk 
luminosity should be only weakly mass de- 
pendent and very nearly the same for the two 
metrics. The binding energy at Vms differs by 
only 3.6% in the two metrics. It is indepen- 
dent of both the mass of the central object 
and r^ns- For the star model used here, disks 
would terminate at r^^ above the surface for 
NSs with masses above 2.1 M© (YM) or 1.7 
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Mq (SM). For stars with sufficiently weak 
magnetic fields (< lO^G) this "gap accre- 
tion" would produce a hard power-law spec- 
trum (Kluzniak & Wilson 1991, Hanawa 1991, 
Walker 1992) arising from bulk flow Comp- 
tonization in a boundary layer. In fact, how- 
ever, radiation (Miller, Lamb & Psaltis 1997) 
and/or magnetic torques (Lai 1998 ) on the 
disk appear to cause the inner disk radius of 
NSs to be larger yet. At present it seems 
likely that kHz QPOs may be generated at 
the sonic point just inside the disk (Miller, 
Lamb & Psaltis 1997). The observed QPO 
frequencies indicate that the sonic radius is 
generally larger than rms, but may approach 
rjns as a limit at sufficiently high accretion 
rates (Zhang, et al. 1998). 

The Keplerian frequency of the innermost 
marginally stable orbit is of considerable in- 
terest. In the YM, Keplerian frequencies are 
given by: 



VGMexp{-2u) 
27rr3/y (1 - u) 



(9) 



For the innermost marginally stable orbit, u 
= 0.191 and r^, = 5.2AGM/c^, this yields 

UK = 2O4OM0/M Hz, compared to 2200 
Mq/M Hz in the SM; a mere 7% difference. 

For neutron stars there is additional energy 
dissipated as accreting particles are brought 
to rest on the surface. Neglecting the ef- 
fects of star spin and evaluating the energy- 
momentum four-vector for zero momentum 
shows the energy of a particle at rest on a star 



of radius R to be E = mQC^^g{R). Thus the 
total accretion energy, external to the star, 
which may be directly radiated to distant ob- 
servers is moc^(l — ^Jg{R)). The fraction ra- 
diated from the star surface itself (neglecting 
fusion, which can be considered as part of the 



binding to the star) would be: 

/, = 1 - - fd (10) 

In the strong- field cases to be considered here, 
over half of the energy release exterior to the 
star occurs within 3R and ~ 90% from within 
lOR. 

(ii) Compact Objects 
A compact object of constant proper density 
(Clapp 1973) is considered here. The details 
are straightforward, but tedious, and have 
been placed in the appendix. The proper den- 
sity is p{r)g{r) = po, a constant. These com- 
pact objects have a characteristic radius and 
mass that are given by: 



V47rGVo 



Mo 



G 



(11) 



(12) 



As shown in Table 3, g(R) reaches 0.25 at the 
surface of a maximum mass NS in the SM. 
Assuming that p{R) = 2.7x10"'^^ g/cm'^, the 
nuclear saturation density, and using g{R) = 
1 - 2u{R) yields po = 1.35x10^^ g/cm^. This 
value of Po yields: ro = 28.2 km and Mq = 
19M0. In turn, these yield a maximum mass 
of 2.45M0 for 14.3 km radius and a maximum 
radius of 15.5 km for the SM. Independent of 
the choice of po, NSs in the upper one third 
of the mass range permitted by the SM would 
have u(R) > 1/6 and have radii small enough 
for gap accretion. 

For the same po in YM, the maximum mass 
is 9 M0 at a radius of 13.5 km. A maximum 
radius of 18.8 km occurs for about 5.1 Mq. 
Mass vs.radius curves are similar in shape to 
those obtained by Friedman & Ipser (1987) 
for their core model. Although the constant 
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Fig. 1. — Radius, R, gravitational mass, m^, 
and surface redshift, z, as functions of free 
baryon mass, m^, contained within a neutron 
star in the exponential metric. Radius is in 
units of tq = cj y^AnGpo and masses in units 
of Mq = c^Tq/G. Heavy line features at the 
lower left are the same quantities for the SM. 



proper density model fails to properly account 
for the softer exterior layers for low mass ob- 
jects, this will be of no concern here. The 
effects of stellar rotation may increase max- 
imum masses (Friedman & Ipser 1987) by 
about 25%, to about 3 M© for the SM and 
11 Mq for the YM. Figure 1 shows compar- 
isons of mass, radius and surface redshift for 
YM and SM. It is interesting that there exists 
a maximum mass for the YM that is depen- 
dent on density. For nuclear density, an 11 
M0 maximum mass may be adequate to cover 
the range of masses that have been found for 
the galactic BHCs. Figure 2 shows properties 
for the SM in more detail. 

Tables 2 and 3 give masses, radii, surface 
potentials, surface redshifts, etc for the two 
metrics. With these quantities, the necessary 
binding fractions of the disk and surface can 
be calculated using equations 7-9. Figure 
3 shows these binding fractions. A profound 
effect of a neutron star surface is implied by 
Fig 3. Up to about ~ 1.4 Mq, the luminosity 
would be fairly evenly divided between disk 
and surface, but for the largest masses in the 
Schwarzschild metric up to 80% of the lumi- 
nosity might arise inside the last stable orbit; 
near or on the star surface. In the YM as 
much as 90% of the luminosity could origi- 
nate from the star surface. The binding frac- 
tions at the star surface can reach 70% for 
large masses in the YM. Large binding frac- 
tions and redshifts would make their surface 
radiations very prominent and very soft. 

(iii) Internal Processes 
After accreted matter reaches the surface there 
is still energy to be produced by fusion and a 
need for additional gravitational binding en- 
ergy to be released for stability of the star. 
In order to properly account for all of the en- 
ergy available, one needs to consider the dif- 
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Fig. 2. — Radius, gravitational mass and sur- 
face redshift as functions of free baryon mass 
for a neutron star in the Schwarzschild met- 
ric. Units are the same as Figure 1. Dashed 
hues show extensions for black holes. 




mg 

Fig. 3. — Efficiencies of conversion of gravi- 
tational potential energy to radiation as func- 
tions of gravitational mass, rrig. Gravitational 
mass is in units of Mq. Solid upper lines are 
(/s + fd)- F^^^ baryon mass, mfe, increases 
from left to right along these curves. The spi- 
ral of the YM occurs because maxima and 
minima of redshift and gravitational mass do 
not correspond to the same free baryon mass. 
Lower dotted line is Z^, which is essentially 
the same in both metrics. Surface radiation 
is represented by the difference between the 
upper and lower lines, excepting black holes. 
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ference between free baryon mass and mass 
bound gravitationally within the star. The 
total gravitational energy release as rest mass 
Amb falls from a large distance and becomes 
bound in the star would be c^(Amb — Amp), 
or c^Am(,(l — drrig/dmb). The term in paren- 
theses being the maximum possible fraction 
of accreted rest mass to be converted to radi- 
ations seen by a distant observer. It is plot- 
ted in Figure 4 as the gravitational binding 
fraction of the object. For a neutron star in 
the SM, it actually diverges shortly after it 
reaches 1.0 just below 2.35 Mq. At this point 
the constant intrinsic density model fails in 
the SM. The addition of more mass simply 
causes a core collapse, though it may be pos- 
sible that some of the binding energy would 
escape through neutrino emission, pair pro- 
duction and gravitational radiation. 

In the YM there is no discontinuity of 
binding fraction but after trapped baryons 
produce a maximum gravitational mass, the 
star can become more tightly bound. The 
more massive neutron stars permitted in ei- 
ther metric swallow baryons but must radi- 
ate their mass equivalents for stability. It is 
apparent that the overall binding fraction is 
much larger than the binding fraction at the 
surface. Some substantial additional radia- 
tions of internal origin should be seen from 
accreting NSs in either metric. It is possi- 
ble that this energy would appear as copi- 
ous photon, pair production or neutrino emis- 
sions, or gravitational radiation. Pair plasma 
emissions might be part of the jets observed 
for some NSs and BHCs (Ghozzi, Bodo, & 
Ghisellini 1998). Whether or not accretion 
energy might be transiently stored has not 
been considered here, but it should be noted 
that mass-energy can be trapped inside the 
photon sphere for u{R) > 1/2 (YM) or 1/3 
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Fig. 4. — Binding energy expressed as a frac- 
tion of the rest mass-energy of accreted parti- 
cles; given by ( 1 - ^^)- Collapse of the core 
in the Schwarzschild metric is shown by the 
diverging line at the left. The curved hue ex- 
tending to the right is for the Yilmaz metric. 
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(SM). Spectral qualities of surface radiations 
should also be affected when the surface lies 
inside the photon sphere, however, this would 
involve analysis well beyond the scope appro- 
priate for this work. 

3. Magnetic Phenomena 

The magnetosphere radius can be estimated 
as the radius for which the impact pressure of 
accreting matter is some fraction, rj, of mag- 
netic pressure (e.g.. Lamb, Pethick & Pines 
1973). For dipole fields and disk accretion, 
with Keplerian motion outside the magne- 
topause, this criterion yields : 
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) 



(13) 



where e is the ratio of radial flow velocity 
to free fall velocity in a spherical flow, 2^ 
is the ratio of disk thickness to radius and 
B is the magnetic field at R. A choice of 
(r^^/e)^/^ = 0.35 brings the magnetospheric 
radius into close agreement with that of the 
much more complicated model of Ghosh & 
Lamb (1992). Using rj^/e = 0.025, M = 
2.SxW^MiAg,B = 10^BgG,R = IQ^Rq cm 
and rh = 10^^mi5 g/s this becomes: 



{78km){ 



Mi.4mi5' 



(14) 



The magnetosphere radius shrinks to accomo- 
date higher impact pressure. It is relatively 
impermeable to diamagnetic plasma, but un- 
stable for curvature perturbations (Rayleigh- 
Taylor). Arons et al. (1984) have described 
shot penetrations of the magnetopause and 
their further fragmentation via Kelvin-Helmholtz 
processes as they fall through the co-rotating 
magnetosphere. At relatively low accretion 
rates individual shots may produce consid- 
erable flickering. The larger surface binding 



fractions for more massive objects in the YM 
might explain the larger amplitude flickering 
observed for BHCs. For recent reviews of 
magnetospheric phenomena, see Campana et 
al. (1998a) and Becker and Triimper (1997). 
It is well-known that pulsars can generate 
very hard spectra via magnetospheric effects 
(Michel 1991, Becker & Triimper 1997) and 
surface impact (Alme & Wilson 1973) as well 
as in shock fronts where magnetospherically 
expelled matter impinges on nebulae. 

The star spin frequency determines the co- 
rotation radius, r^, for which the magneto- 
spheric equator rotates at the Keplerian or- 
bit frequency. In (SM), = (G'M/47rV|)^/l 
When the magnetospheric radius hes outside 
the co-rotation radius, accreting matter re- 
ceives a super-Keplerian push at the magne- 
topause. It is swept outward and is unable to 
penetrate the magnetosphere until a bound- 
ary layer with sufficient pressure builds. This 
is the well-known "propeller effect" (lUari- 
anov & Sunyaev 1975). When the magneto- 
sphere is forced inside the co-rotation radius, 
the spectrum changes from a hard, low state 
to a state of higher luminosity and softer spec- 
trum due to greater surface contributions. If 
the transition radius is large there may be 
a large flare as mass strikes the star. Since 
the co-rotation radius is determined by the 
spin, there is a strong luminosity-spin corre- 
lation for these ffares (Stella, White & Rosner 
1986). Additional evidence has recently been 
found (Cui 1997) for propeller effects in the 
accreting x-ray pulsars, GX 1-1-4 and GRO 
J1744-28. Sharp transitions were observed at 
luminosities of ~ 10^^ erg/s, for which x-ray 
pulsations ceased, but became re-established 
at slightly higher luminosity levels. Zhang, 
Yu & Zhang (1998) and Campana, et al. 
(Campana, S. et al. 1998b) have indepen- 
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dently analyzed the decline of the 1997 out- 
burst of Aql X-1 and found strong evidence of 
a propeller effect and an accompanying spec- 
tral state transition with very little luminos- 
ity change. The transition was followed by 
a rapid cutoff of luminosity that was arrested 
only briefly as radiation pressure from the sur- 
face emissions ceased. 

For sufficiently large rc, the maximum lu- 
minosity that can be achieved without sur- 
face contributions corresponds to ~ Tc and 
Lc = GMm/2rc: 



eGM 



(15) 



or 



Le = {1.2xlO^'^erg/s)Bg^Re^i^2^MiA~^ (16) 

where us = lOOz/2 Hz. Campana et al. (1998) 
have produced the same equation for spher- 
ical accretion except for a multiplier of 20 
in place of the 1.2 used here. The smaller 
multiplier merely reflects the smaller magne- 
tospheric area impacted by a disk. 

For objects more massive than ~ '^Mq, 
classical physics expressions for surface lumi- 
nosity are inadequate, however, we can use 
L — (yfs + fd)i^c^i where y is the fraction of 
mass reaching the magnetopause that pene- 
trates to the surface. Using this to evaluate 
fn and setting ~ in Eq. 13 provides a 
relation of the luminosity for accretion at the 
co-rotation radius to the magnetic fleld and 
spin of the star: 



V2e{GM)^ 



(17) 



or 



L = {9xl0^'erg/s){yfs + U)BlRlviMA 



2 d6,.3 H/f 3 



Eqs. 16 and 18 should yield identical results 

for y = 0, however, since the classical and 
relativistic binding energies differ slightly for 
the inner disk, there can be spurious small 
differences. 

Lmin will be used to designate the accretion 
rate for which y=l, r = and all of the 
accreting matter reaches the star. If the co- 
rotation radius is relatively large compared to 
the star radius, large luminosity changes can 
occur at the transition from Lc to I/mm, even 
if the accretion rate changes very little. Even 
slight departures from y = might produce 
such strong radiation pressure that the full 
transition to Lmin is prevented. 

The magnetosphere is limited by the light 
cylinder radius, r^, at which a co-rotating 
magnetosphere equator would move at the 
speed of light. For fast spinners with us ~ 
400i7^, Tic can be less than ~ 120/cm. When 
the magnetospheric radius lies beyond r/c the 
radio pulsar mechanism may become active. 
Both pulsed point source x-ray emissions and 
diffuse harder radiation at lower luminosity 
from a synchrotron nebula have been ob- 
served for the isolated 7-ray pulsar PSR 1055- 
52 (Shibata et al.l997 ApJ 483, 843). 

The maximum quiescent luminosity that 
can be produced solely by disk accretion cor- 
responds to Tm — Tic- With L — GMm/2r, 
this yields: 



^r](27rus)^VGMB^R^ 



■'q,max 



2V2i 



ec2 



(19) 



or 



{3AxlO^"erg/s)Mi,jBg'Re''u2 



^q,max 



(20) 

If a pulsar wind drives the quiescent luminos- 
ity, the rotational energy loss rate is (Bhat- 



(18) 
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tacharya & Srinivasan 1995): 

E = A.'kHvsvs = ^ (21) 

where I is the moment of inertia of the star. 
The right member of this equation is empiri- 
cal, lacking magnetic axis inclination depen- 
dence and perhaps factors of ?7 and e, however, 
it and the middle member have been taken as 
an operational definition of the magnetic mo- 
ment, BB?, for NSs with I = W'^^g - cw?. 
Becker and Triimper (1997) have provided 
a tight correlation between x-ray luminosity 
and rotational energy loss rate for which the 
x-ray luminosity is 1.0xlO~^47r^/^'5i/5 again 
with fixed I = 10^^ g cm?. While most of 
the correlated sample likely consists of stars 
near IAMq, what has really been correlated 
is x-ray luminosity and ps^s- Had a different 
value of I been used in the calculation of E, 
the multiplier that would correlate with the 
x-ray luminosity would be 1.0xlO~^a;10*^^//. 
For use with the YM, / = 2.85a;10^^ g - cm^, 
appropriate for IAMq and Rq — 1.5 will be 
used here. In this way the YM model can be 
satisfied while maintaining relatively accurate 
x-ray luminosities. Combining this with the 
right member of Eq. 21 yields: 

Lq = (1.3x10^° erg/s)BlRlv^ (22) 

In contrast to previous equations, this is merely 
a well grounded empricicism that represents 
some loosely defined x-ray luminosity rather 
than a bolometric luminosity. Still it is to be 
preferred over Lg^max for quiescent luminosity. 

For accretion to reach the star surface, 
with quiescent luminosities of 10^^ — 10^'^ erg/s 
being above Lc, would require a combination 
of very weak field and small co-rotation ra- 
dius (fast spin). If most NSs reach a spin 
equilibrium for long time average luminosi- 
ties near those in accord with Eq. 18 with y 



= 1, (White & Zhang 1997) then their quies- 
cent luminosities are highly unlikely to derive 
from either surface accretion or accretion at 
the light cylinder radius. 

Becker and Triimper 's (1997) strong corre- 
lation of x-ray luminosity with magnetic field 
strength and spin rate down to luminosities 
well below 10^° erg/s strongly suggests that 
a preponderance of the quiescent luminosity 
is magnetospheric. There are reasons for be- 
lieving that the power-law portions of qui- 
escent emissions originate near the magne- 
topause (see below). Nevertheless, some qui- 
escent luminosity may derive from a cooling 
NS surface. Excess soft radiation correspond- 
ing to kTjji, Ki 0.1 — 0.2 keV and correspond- 
ing to small ~ 10 — 100/cm^ areas has been 
found to be a part of the quiescent emissions 
of both NSs and BHCs (Tanaka & Shibazaki 
1996). It is conceivable that these could be 
polar cap emissions of NSs, but Heindl & 
Smith (1998) have shown that proper atten- 
tion to flourescence and reflection substan- 
tially reduces the apparent temperature of the 
excess soft emissions. Further, even the re- 
maining apparent temperature may be larger 
than the surface temperature of a NS, de- 
pending on its atmosphere. As shown in dips, 
some of the soft excess also appears to orig- 
inate from an extended disk as well. Thus 
the origins and magnitudes of contributions 
to the quiescent luminosity are somewhat un- 
clear at the present time. For this work, 
for luminosities generally above lO'^^ erg/s, it 
will be assumed that the power-law magne- 
tospheric emissions dominate the hard states. 
Ultrasoft blackbody {kTbi, = 0.048 keV) sur- 
face radiation has been observed for the iso- 
lated gamma-ray pulsar, Geminga (Halpern 
& Wang 1997). No evidence of polar cap heat- 
ing was found. Non-thermal power-law mag- 
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netospheric emissions in the spectral range 0.5 
to 2 keV were present. Obviously, ultrasoft 
and power-law spectral components are not 
signatures of black holes. Geminga is a high 
energy gamma-ray source with a l.dxlO^^G 
magnetic field. It is the only known radio- 
quiet, rotation-powered pulsar with a strong 
magnetic field, however, some of the BHCs 
might also have strong fields; see below. 

Brown, Bildsten & Rutledge (1998) have 
found that crustal heating during accretion 
episodes resets the core and surface tempera- 
tures of NSs and might allow them to produce 
quiescent luminosities of as much as ~ 
erg/s. An actual temperature of as little as 
~ 0.03 keV, would contribute ~ 10^^ erg/s to 
the quiescent luminosity. This would produce 
a radiant equilibrium temperature of about 
5000 K for optically thick matter at a distance 
of about 10^ km. The inner accretion disk 
would need to be beyond this for thermal- 
viscous stability. For spins in excess of 0.005 
Hz, the light cylinder radius is smaller than 
10^ km. This effectively eliminates Lg max for 
the production of quiescent luminosity. It also 
sets a minimum inner radius at the start of 
outbursts. The inner disk must fill on a vis- 
cous timescale, producing a delay of days be- 
tween the onset of optical and x-ray luminosi- 
ties. Thus the 6 day delay (Orosz, et al. 1997) 
between the onset of optical and x-ray lumi- 
nosity increases for the 1996 outburst of GRO 
J1655-40 would be expected. Hameury, et al 
(1997) have calculated 6 days as the interval 
for viscous infiow from 2a; 10^ /cm. Delays for 
refill of the disk have been observed for dwarf 
novae and for the NS, 4U 1608-52. 



4. The Magnetic Fields of NSs and 
BHCs 

As additional information is acquired, some 
of the tentative conclusions drawn here will 
surely need alteration. The examples chosen 
are those for which I have sufficient informa- 
tion to initiate a discussion. They do not all 
fit neatly into any current scheme. 

(i) Neutron Stars 
It is well established that radio pulsars may 
have magnetic fields ranging from 10^ G for 
millisecond pulsars to 10^'^ G for slower spins. 
If not for spins revealed by bursts or kHz 
QPO frequency differences, it would be diffi- 
cult to determine the magnetic fields of Atolls 
and Zs. They are not pulsars in any normal 
sense. Nevertheless, ratios of luminosities can 
be used to determine their spins and magnetic 
moments. Since all of the luminosity equa- 
tions considered here depend on the dipolc 
moment, BR\ it can be eliminated by the 
use of luminosity ratios, thus permitting the 
spin to be estimated. Once the spin is known, 
B can be calculated. For consistency with the 
constant proper density star model, R = 15 
km for a 1.4 Mq will be used here. Since the 
quantities being determined from luminosity 
ratios are actually spins and dipole moments, 
an overestimate of R may cause an underes- 
timate of B without other effect. In order to 
use the ratio method, two of L^^^, L^, Lq^^^ax 
or Lq must be identified. 

Aql X-1 is an important test case for the 
use of ratio methods. Campana et al (1998b, 
Fig. 1) reported spectral hardening beginning 
at Lmin = 1.2a;10^^ erg/s and complete cessa- 
tion of the rapid decline at about 1.2x10^^ 
erg/s. Identifying these as Lmin and Lq^max, 
the ratio yields ug = 512 Hz for mi.4 = 
1. This agrees moderately well with a spin 
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frequency of 549 Hz observed during burst 
QPOs. Using the QPO spin frequency and 
Eq. 18, a magnetic field of 1.3x10^ G can 
be calculated. Aql X-1 displays ultrasoft and 
power-law components in its quiescent spec- 
trum with the power-law component account- 
ing for more than half the luminosity. Using 
549 Hz and LSajlO*^ G, the magnetospheric 
component of quiescent luminosity is calcu- 
lated to be 2.3x10^^ erg/s. Garcia, McChn- 
tock & Narayan (1997) report a quiescent lu- 
minosity of 4a;10^^ erg/s. Further, using 549 
Hz and 1.3x10^ G permits the calculation of 
Lc = 3.8x10^^ erg/s, as obtained by Gampana 
et al (1998). 

SAX J1808. 4-3658 A second test case is 
provided by the recent outburst of the 401 Hz 
x-ray pulsar and burst source, SAX J1808.4- 
3658, (Gilfanov et al. 1998, Hcindl & Smith 
1998) which displayed a light curve very sim- 
ilar to Aql X-1. It reached a luminosity 
level of 5x10^^ erg/s, dechned slowly to about 
1.5x10'^'^ erg/s and then declined rapidly with 
a brief arrest at 4x10^^ erg/s before reaching 
about 10^^ erg/s. Between the start of rapid 
decline and the brief arrest, 401 Hz pulsations 
changed from 4% rms to undetectable, thus 
indicating that surface impacts ceased. Iden- 
tifying the luminosities at the start of rapid 
decline and brief arrest as Lmin and allows 
the calculation of a magnetic field strength of 
2.1x10^ G, a convincing spin frequency of 419 
Hz, Lg^rnax = 8x10^^ erg/s, and Lg = 1.7x10^^ 
erg/s. 

When the accretion rate is declining, the 
onset of the propeller effect may be followed 
by a rapid cutoff of luminosity. Although 
many NSs (and BHGs) have shown a steep 
decline of luminosity from the vicinity of 10^^ 
erg/s, this may not always be clear evidence 
of the onset of a propeller effect. Without 



the cessation of pulses it would have been un- 
convincing to identify the propeller effect for 
SAX J1808.4-3658. Observations of Type I 
bursts at luminosities above 10^^ erg/s show 
that accreting mass reached the surface. Yet 
the spectrum remained a hard power-law ex- 
tending to 100 keV over two decades of lu- 
minosity change. A detailed spectral analy- 
sis (Heindl & Smith 1998) revealed the pres- 
ence of iron fiourescence and a small, not very 
well constrained, soft excess which was mod- 
eled MBB. An additional unusual feature 
of this source is the soft x-ray lags reported 
by Gui, Morgan and Titarchuk (1998). An 
obvious possibility for explaining the unusual 
spectral and timing characteristics would be a 
very small inclination view with strong polar 
axial magnetospheric emissions obscuring the 
surface. Alternatively, there might be some 
unusual surface mechanism that produced the 
hard spectrum. At any rate, the presence of 
only hard spectral components cannot be in- 
terpreted lack of surface accretion. 

The spin frequencies determined from lu- 
minosity measurements for Aql X-1 and SAX 
J1808.4-3658 are in convincing agreement with 
burst QPO frequencies, however, if their as- 
sumed masses had been doubled, their cal- 
culated spins would have been halved and 
their magnetic fields roughly doubled. For 
4U 1636-54 and KS 1731-26, the frequency 
difference of kHz QPOs has been near half 
of the burst QPO frequency (Swank 1998). 
The question arises as to what is the true 
spin frequency. In this case it seems to be 
the burst QPO frequency. It has been sug- 
gested that the burst QPO frequency might 
be double the actual spin frequency due to 
two antipodal hot spots on the star. It seems 
improbable that bursts would be simultane- 
ously initiated at two different places on the 
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star surface (Strohmayer et al. 1998). On 
the other hand, the detonation of some 10^^ 
megatons at one pole might quickly ignite an- 
other pole less than 50 km distant, depending 
on the shielding afforded by an intervening 
stellar mass (Miller 1998). 

4111608-52 provides another test case for 
the use of ratios. Spectral hardening be- 
gan in decline at Lmin = 10'^'' erg/s (Mit- 
suda et al. 1989). Garcia, McClintock & 
Narayan have reported a quiescent luminos- 
ity of Lg = 2x10^^ erg/s, while Tanaka & 
Shibazaki (1996) report 6x10^^ erg/s. For 
these luminosities, the ratios L„iin/Lg and 
Eqs. 18 and 22 yield spins of 566 Hz and 
275 Hz. The frequency difference of twin kHz 
QPOs for this source is near 230 Hz, however, 
this might be only half of a 460 Hz actual spin 
frequency. The conflicting luminosity data 
leave this case unresolved, but it illustrates 
an important point; that the results are prob- 
ably correct to within small numerical factors 
anyway. The magnetic fields estimated from 
Ljnin and the two QPO spins are 4.6x10^ G 
and 8.4x10^ G. The higher spin and weaker 
field are the more hkely attributes of an Atoll 
source. 

4U 1820-30 apparently has a mass in ex- 
cess of 2 M0 based on its kHz QPO limit 
frequency of 1060 Hz (Zhang, et al. 1998). 
Based on QPO frequency differences it spins 
at either 275 or 550 Hz. Both of the high fre- 
quency QPOs were visible in the low state 
at a luminosity of 1.5x10^'' erg/s but cut 
off when the spectrum softened for luminosi- 
ties in excess of S.lxlO^^ erg/s. Assuming 
that the start of twin QPOs would corre- 
spond to mass reaching the surface implies 
Lmin = l.SxlO^^ For 2 Mq, Eq 18 yields ei- 
ther 5a;10*^ G or 1.1x10^ G. The higher field 
would be an unusually strong field for an 



Atoll source, but perhaps not for one of 2 
Mq. Calculated quiescent luminosities for the 
two cases are 5x10'^'^ erg/s and 1.5x10^^ erg/s. 
Spin and magnetic field estimates for other 
Atoll sources, 4U 2129 + 47, KS 1731-26 and 
4U 1730 - 335 are included in Table 1. 

Cir X-1: Neither spin rate nor mass are 
available for the enigmatic Cir X-1. Lmin = 
6.3x10^^ erg/s and = 3x10^^ erg/s can be 
determined from observations just before and 
two hours after a spectral hardening transi- 
tion on Sept 20-21, 1977 (Dower, Bradt & E 
Morgan 1982). The ratio permits an estimate 
of a 32 Hz spin followed by a magnetic field 
estimate of 8x10^° G. These results were ob- 
tained for M1.4 = 1.0 and Rq — 1.5. If the 
mass of Cir X-1 were IMq, the spin would be 
38 Hz and the magnetic field 8.2x10^° G. The 
luminosity ratio method is relatively insensi- 
tive to mass, provided that the correct accre- 
tion efficiency is used. The calculated quies- 
cent luminosity from Eq. 22 is (0.9 — 6)a:;10^^ 
erg/s, which is well below detector sensitivi- 
ties. Cir X-1 experiences a much wider range 
of accretion rates on orbital time scales than 
most NSs. Bradt, Shirey and Levine (1998) 
have provided a fascinating series of hardness- 
intensity diagrams for this source. 

The Z and Atoll NSs of Table 1 have rela- 
tively small differences between L„,j,„ and 
due to their rapid spins. The ratio of these, 
A, should be: 



yfs + fd 

fd 



(23) 



To first order with y = 1, r = r^s in the 
numerator and r = Tc in the denominator, 
A fti 2u{R)/u{rc) ^ 2rjR a; 2 3 for 
fast spins. This result has been discussed for 
spherical accretion by Corbett (1996). The 
bolometric luminosity change when the mag- 
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netopause reaches the co-rotation radius may 
be much smaller than implied by the ratio 
2rjR iiy <l. 

If NSs of the fast spin Atoll and Z classes 
have been spun up to high equilibrium spin 
rates by accretion, then the long-term average 
luminosity would be near Lmin- White and 
Zhang (1997) used these luminosities and the 
known spins of several Atoll and Z class NSs 
to estimate their magnetic fields. Values of 
Lmin or Lc found for the spectral state tran- 
sitions of several other NSs have been used 
here to determine some of the NS magnetic 
fields shown in Table 1. The results, though 
intentionally sparse, fit the pattern suggested 
by van der Klis (1994) and found by White 
and Zhang (1997). Z sources appear to have 
fields that are stronger than those of Atolls. 
There appears to be a Lmin ol correlation 
in the results. White and Zhang suggested 
that the correlation might be spurious; that 
there might be comparable magnetic fields in 
all cases with a magnctospheric radius that is 
more weakly dependent on accretion rate for 
Zs where radiation pressure dominates the in- 
ner disk. The clear differences of properties 
of Atolls, Zs and Cir X-1 suggests that the 
correlation represents a progresssion of a real 
physical property. 

Kilohertz QPOs have now been observed 
for many NSs, including both Atoll and Z 
types. The Atoll source 4U 1820-30 showed a 
limit frequency of 1060 Hz, which appears to 
correspond to the innermost marginally sta- 
ble orbit for ~ 2Mq (Zhang, et al. 1998). 
This mass estimate obtained from the limit 
frequency of the QPO must be taken seri- 
ously. Cyg X-2, a Z source with a mass > 
I.9M0 found from orbit parameters (Casarcs, 
Charles & Kuulkcrs 1998) has produced kHz 
QPOs of 1007 Hz but without reaching a lim- 



iting frequency; perhaps implying a mass near 
2M0. The BHC GRO J1655-40, at 7 Mq 
(Bailyn et al. 1995) has produced a 298 Hz 
QPO while in a high luminosity state (Remil- 
lard et al 1997). If this is the Keplerian fre- 
quency of the marginally stable orbit it im- 
plies a mass of 7.4 M© (SM) or 6.8 M© (YM). 
One clear implication of these results is that 
the accretion disks of both NSs and BHCs 
may reach the innermost marginally stable or- 
bit. With binding fractions hmited to about 
6% in all such cases, the differences between 
NSs and black holes must be found inside the 
last stable orbit, in the differences between 
surfaces and event horizons. The disks, ex- 
cepting differences due to irradiation from the 
central source, should be much the same, (ex- 
cepting those of Kerr metric black holes). 

If kHz QPOs arise from clumps of mass 
falling from the inner disk radius to the sur- 
face, then the disk inner radius must be less 
than the co-rotation radius. Radii of Keple- 
rian orbits at the frequencies at which kHz 
QPOs first appear are generally about two 
thirds of the magnctospheric radii determined 
from Eq. 13. Dong Lai (1998) has sug- 
gested that the sonic radius might be identi- 
cal with the magnetospherc radius when out- 
side the innermost marginally stable orbit. 
This might be accommodated here by suit- 
able choice of the impact pressure-magnetic 
pressure balance. To do so would reduce the 
magnetospherc radii calculated from Eq. 14 
and change the calculated dipole moments de- 
termined from the luminosity equations with- 
out other effect. 

(ii) Black Hole Candidates 
If BHCs are really massive NSs that show 
neither pulses nor bursts then spectral state 
switches and ratio methods must be used to 
determine their spins and magnetic fields. As 
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is the case for Atolls and Zs, coincident spin 
and magnetic axes arc necessary to account 
for the absence of x-ray and radio pulses. For 
BHCs, it is possible in some cases to treat 
the inner disk radius determined from MBB 
fits to spectra as a co-rotation radius. If the 
massive NS hypothesis is correct, then the 
substantial surface contributions may explain 
the small inner radii that are often found for 
MBB fits to BHC high state spectra. For low 
states, however, in which surface accretion is 
largely blocked by propeller effects, the BHC 
spectral models are quite applicable to NSs. 
Table 1 lists estimated spin and magnetic field 
parameters for the BHCs discussed below. 

GRS 1124.-68 shows characteristics that 
clearly illustrate the nature of the spectral 
state switch. Misra & Melia (1997) showed 
that its inner disk radius increased to about 
27Rg {Rg = 2GM/c^, 400 km for 5 Mq) af- 
ter the spectral state transition. Ebisawa, et 
al (1994) give a luminosity of Lc = 6.6x10^^ 
erg/s for the low state. The assumption of 
Tc = 400 km permits estimates of Us = 15 Hz 
and B = l.lxW^ G. In turn L^i^ = 2AxlO^^ 
erg/s. Eq 22 yields Lq = 3.9x10^^ erg/s. 
The calculated L^j„ fortuitously agrees ex- 
actly with observations. The observational 
upper limit for quiescent luminosity is 4x10^^ 
erg/s, which is consistent with Lg calculated 
for the pulsar wind mechanism. 

Zycki, Done and Smith (1998) fitted the 
spectrum of GRS 1124-68 with a power- law 
plus a reprocessed fraction that included an 
iron line, and a relatively cool disk extend- 
ing from Tin to bxlO'^Rg (with Rg as de- 
fined above). A soft thermal component was 
modeled separately by optically thick Comp- 
tonization of soft seed photons from a MBB. 
By separately accounting for soft components 
their model would not have r,„ forced to a 



small value corresponding to a star surface, if 
one existed. They found little, if any, change 
of Tin at the transition to the intermediate 
spectral state. There was a dramatic decrease 
of the degree of ionization of the disk which 
they attributed to the disappearance of the 
soft component. In addition, they attributed 
a hardening of the power-law component to 
the decreased availability of soft photons for 
Compton coohng. A strong reflected compo- 
nent was found all the while the low state de- 
cline continued, which indicates that optically 
thick material was present within about 700 
km throughout the observations that were an- 
alyzed. In all respects, these results exactly 
accord with expectations for a massive NS. 

GS 2000 + 25: Zycki, Done and Smith 
(1997a) have found a similar spectral state 
switch for GS 2000 + 25. When a strong soft 
component is present the disk is highly ion- 
ized and iron flourescence strongly relativisti- 
cally smeared. The disappearance of the soft 
component is accompanied by hardening of 
the power-law and a decrease of ionization. 
Little change of inner disk radius occurs at 
the transition, but the subsequent decline of 
luminosity is accompanied by an increase of 
inner disk radius. Using an inner radius of 
(50 - 100) Rg (1040 - 2080 km with Rg = 
2GM/c^) as the co- rotation radius (Done, C. 
personal communication) and Lc — 1.5x10^^ 
erg/s for the hard state of Dec. 16, 1989, 
Us = (4.5 - 1.6) Hz and B = (1.4 - 6.5)xl0^^ 
G are obtained. The calculated quiescent lu- 
minosity is (3.7 — 1.3)a;10^° erg/s; well below 
the observational upper limit of 2x10^^ erg/s. 

GS 2023 + 338 never displayed an ultra- 
soft spectral component, but on May 30, 1989 

it changed luminosity by a factor of 21 in 10 s, 
from a saturated 10'^^ erg/s to 4.8x10"^^ erg/s 
(Tanaka & Lewin 1995). 10^^ erg/s is the 



17 



Eddington limit for 7 Mq. The iron floures- 
cence line was not obviously present in the 
high state but showed up prominently at the 
lower luminosity. This might be caused by 
a distended inner disk covering the magneto- 
sphere. When the luminosity diminished, the 
reduction in the 1-10 keV band was greater 
than for higher energies. These phenomena 
strongly suggest that the magnetosphere ex- 
panded beyond the co-rotation radius when 
the luminosity dropped. Assuming this to be 
the case, apparently A a; 21. For 7Mq and 
A = 21, u{R) = 0.56, /, + fdirras) = 0.428 
(Table 2), and /d(r) = 0.02. Solving for u(r) 
yields 0.0434, which then gives r = 240 km 
as the CO- rotation radius. Using m = L/ f^c^ 
with L — 4.8x10^^ erg/s gives m — 2.7x10^^ 
g/s. Solving Eq. 14 for the magnetic field 
then yields 9. 5a; 10^° G. The spin frequency 
(YM) determined from the co-rotation radius 
is 39 Hz. 

From spectral analysis Zydd, Done and 
Smith (1997b) found rj„ = 263 km and 368 
km, respectively, for June 20, 1989 and July 
19-20, 1989. Luminosities (0.1 - 300 keV) 
for these dates based on spectra of Tanaka 
(1992) are 1.9a;10=*^ erg/s and 6.3x10=^^ erg/s, 
respectively. These radii give values of fa of 
0.018 and 0.0134 and values of m are then 
1.2a;10i^ g/s and 5.2x10^^ g/s. Values of B = 
7.3 xlO^° G and 8.8x10^° G then follow from 
Eq. 14. The three independent estimates of 
magnetic field are about as consistent as can 
be expected with an Eddington limit affect- 
ing the first result. Stated another way, if the 
9.5x10^° G field is used with the values of m 
to calculate the magnetospheric radii, values 
of 286 km and 365 km are found. These are 
within 9% of spectral fit values that have sta- 
tistical uncertainties in excess of 50%. 

Using the 39 Hz spin, the quiescent lu- 



minosity would be (5.2 — 8.9)xl0^^ erg/s for 
(7.3 — 9.5)a;10^°G. Quiescent luminosities of 
6.9x10^'^ erg/s Chen, Shrader & Livio (1997), 
8a;10^^ erg/s (Tanaka & Shibazaki 1996), and 
1.6x10^^ erg/s (Garcia, McChntock and Narayan 
1997) have been reported. In any event, the 
calculated quiescent luminosity depends very 
strongly on spin, which is quite uncertain in 
this case with Lmin so near the Eddington 
limit. 

IEI74O.7-2942 is a jet source. Vilhu et 
al. (1997) fitted a "sombrero" model to a low 
state spectrum and found r.j„ ~ 7 — 8Rg for 
the IOMq they assumed. Taking this to be 
an inner radius of Vc = 225 km and to 
be the corresponding 3x10^''' erg/s for 7Mq 
yields us > 43 Hz and B > 6.9x10^° G. The 
hard spectrum with a small inner radius is in- 
compatible with the advective accretion flow 
(ADAF) model of Narayan, Garcia and Mc- 
Chntock (1997). The calculated quiescent lu- 
minosity is < 7x10^^ erg/s. 

Cygnus X-1 has been observed to emit 
gamma rays (Ling et al. 1987). It ex- 
hibits spectral state switches with very lit- 
tle change of bolometric luminosity at about 
5x10^^ erg/s (Belloni et al. 1996). Misra & 
Melia (1997) obtained an inner disk radius of 
56.7 Rg from spectral fltting of an interme- 
diate (hard) state. For 10 Mq, this is 1700 
km. Assuming this to be the co-rotation ra- 
dius yields a spin frequency of 2.6 Hz. Using 
R = 14 km in Eq 18 then yields B = 1.5x10^^ 
G. Magnetic flelds of this magnitude have 
been associated with other gamma emitters 
such as the Crab and Vela pulsars. A lack of 
pulses would require nearly coincident rota- 
tion and magnetic axes and a lack of multi- 
pole moments. These are stringent, but not 
impossible, requirements. The calculated qui- 
escent luminosity, should it ever be observed 
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for comparison, is 1x10^^ erg/s. 

A0620 00: Spectral hardening began about 
100 days after the start of the 1975 outburst 
and continued until interrupted by a reflare. 
Lmin = 8xW^^ erg/s (E. Kuulkers 1998) and 
Lg = 10^^ erg/s (Garcia, Narayan & McCHn- 
tock 1997) along with Eqs 18 and 22 permit 
the estimates, us = 13.5 Hz and B = 2.3x10^° 
G. The kinetic energy that the magnetosphere 
can impart to impinging matter is generally 
much less than the binding energy for the in- 
ner disk, but not for the light cylinder. This 
may be involved in the reflare process. In the 
case of Atolls and Zs, the spins are so fast that 
matter is too tightly bound at the light cylin- 
der, hence little reflare is observed for them. 

GRO J1655-40: L^in was reached at 
3.1a;10^''' erg/s approximately July 29, 1996 
(Mendez, Belloni and van der Klis 1998). The 
rapid decline was arrested at 1.5x10^^ erg/s, 
which can be taken to be Lc- Eqs. 16 and 
18 then yield ug = 34 Hz and B = 2.1x10^° 
G. The calculated quiescent luminosity for 
these values exactly matches the observed 
2.5x10^^ erg/s reported by Garcia, McClin- 
tock & Narayan (1997). 

GRO JO422 + 32 never showed an ul- 
trasoft peak and stayed in a spectrally hard, 
gamma emitting state during its 1992 out- 
burst (Grove, et al. 1998). Lacking evidence 
of any spectral state change, only possible 
limits for spin and field can be estimated here. 
Using 3.6M0, (Chen,Shrader & Livio 1997) 
maximum luminosity 7.9x10^^ erg/s as Lmin 
and quiescent luminosity < 7.9x10^^ erg/s 
(Garcia, McClintock and Narayan 1997) as 
limits, one obtains vs < 0.36 Hz and B > 
1x10^^ G. Such large magnetic fields may oc- 
cur in nature, but it is possible that the spec- 
trum might have remained hard for the same 
reasons, whatever they were, as SAX J1808.4- 



3658. A slow spin would be consistent with 

the remarkable lack of high frequency features 
in the power density spectrum. This lack is 
at odds with accretion models that produce 
all hard x-rays and gamma emissions from 
regions near the innermost marginally sta- 
ble orbit (Grove et al. 1998). Grove et al. 
found a strong QPO, independent of photon 
energy or luminosity at 0.23 Hz. If this is the 
spin frequency, the magnetic field could be as 
high as 2x10^^ G. Another possibihty is that 
the rapid decay after about 10^^ erg/s might 
correspond to the onset of the propeller ef- 
fect. This would still imply a large magnetic 
field of 10^^ G. In either of these slow spin 
cases, the quiescent luminosity would be be- 
low 5x10^^ erg/s; still below the observational 
upper limit. 

GRS 1915 + 105 displays oscillations with 
peaks above the Eddington limit followed by 
hard states that are lower in luminosity by a 
factor of 3 - 6. The conventional view is that 
the inner disk becomes unstable and a flare 
is produced by unspecified means as it falls 
into a black hole. Belloni et al. (1997) have 
attributed the intervals between fiares to the 
time required for the inner disk to refill on 
a viscous time scale, which is likely correct. 
There is, however, a more credible explana- 
tion of the flares. They could arise from a 
star surface. If plasma falling from the in- 
nermost marginally stable orbit can produce 
a flare on the surface that is beyond the Ed- 
dington limit, radiation pressure could then 
push the disk back beyond this orbit and tem- 
porarily cut off the flow to the surface. After 
the inner disk refills, another fiare follows, etc. 
MBB spectral fits (Belloni et al 1997) show 
that Tin oscillates between about 20 km and 
80 km, but occasionally reaches only 55 km, 
followed immediately by another burst. 55 
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km would be the innermost marginally stable 
orbit radius (YM) for TMq, which is therefore 
adopted as the mass. 

A 67 Hz QPO has been observed (Remil- 
lard et al 1997, Remillard & Morgan 1998) 
that is sharp [Q > 20), stable for factors of 5 
luminosity change and stable over six months 
time. It is observed most clearly for lumi- 
nosities for which surface radiations should be 
seen, if there is a surface. Speculating that 
this is the spin frequency, a co-rotation ra- 
dius of 163 km (YM) is imphed. The specu- 
lative 67 Hz spin frequency must be tested by 
further observations. There are two sources 
of drifts and broadening for this QPO. Or- 
bital characteristics may imprint cyclical vari- 
ations on it. Observing it at high luminosi- 
ties with radial expansion of the surface layers 
may broaden it. If some of the hard secondary 
bursts that have been observed (Taam, Chen 
and Swank 1997) originate on the surface, as 
seems to be implied by the small radii deter- 
mined for them, they may also show identifi- 
able spin QPOs. 

The low state between flares reached about 
4a;10^^ erg/s for a radius of about 80 km. 
Evaluating for 80 and 163 km, calculat- 
ing m = 9.3x10^^ g/s for the 80 km radius 

-7 

and scaling it for 163 km, (a Vm ) permits 
the calculation of = 1.9x10^''' erg/s. If the 
accretion rate of 9.3,t10^^ g/s for the cutoff 
state were to impact the surface after reflUing 
the inner disk, the luminosity could be driven 
to 3.6x10^^ erg/s, if not for the Eddington 
limit (~ 10^^ erg/s). The excess energy above 
the observed 1.6x10'^^ erg/s may well be what 
drives the jets. In the 10 - 50 s durations of 
flares, excess energy of lO^'^ — lO^-*^ erg might 
be imparted to jets. The magnetic held, esti- 
mated from Lc and Eq. 16 is 2.8x10^'^ G. It 
is noteworthy that this exceeds the minimum 



10^ G found by Gliozzi, Bodo & GhiseUini 
(1998) to be necessary for the fleld strength 
at the base of the jets of GRS 1915+105. 

Eq 22 yields an expected qTiicsccnt x-ray 
luminosity of 7a;10^^ erg/s. In all parameters, 
GRS 1915 + 105 seems to be like GRO J1655- 
40 and Cir X-1, both of which are also chaoti- 
cally variable in both x-ray and radio regimes. 
It should be clear from the low value of L^, 
that accreting mass should have reached the 
surface at all times for the observed oscilla- 
tion luminosities under consideration. That 
it apparently did not, judging by the hard 
spectrum after a burst, is due to the empty 
inner disk after the burst. The star behaves 
as a relaxation oscillator between the super- 
Eddington flares and innermost marginally 
stable orbit. For oscillations, it is critical that 
the magnetic fleld be of the right strength to 
just require flows near the Eddington limit to 
drive the magnetosphere near the marginally 
stable orbit. Massive stars with very large 
surface efficiencies can then provide a strong 
burst as mass hits the surface. The oscillating 
state might require a balance of fields and ac- 
cretion rates too delicate for achievement by 
low mass neutron stars that have no accretion 
gaps. 

5. Spectra of Massive NSs 

For this discussion it must be remembered 
that there are only two reliable distinguish- 
ing spectral characteristics of BHCs. These 
are the simultaneous presence of hard spectral 
tails and ultrasoft spectral components above 
20^7 erg/s and larger amplitude flickering for 
some BHCs. Both are indicators of larger 
masses for BHCs in the massive NS model. 
No event horizon explanation has been pro- 
posed for them for black hole models. Flick- 
ering can be explained by shot penetrations 
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of the magnetopause, as illustrated by the 
flickering pulsar, V0332+53. A shot mech- 
anism has the virtue of being able to account 
for the strong coherence (Nowak et al. 1998) 
of soft and hard emissions during flickering. 
Stronger flickering of the BHCs is due to their 
larger surface binding energies. The quench- 
ing of the hard spectral tail of the low mass 
NSs is due to the fact that 10'^'' erg/s is a 
much larger fraction of the Eddington limit 
for low mass NSs. At this and higher lumi- 
nosities, radiation pressure inflates the inner 
disk and the flood of soft photons cools the 
Comptonization region. The magnetosphere 
would not be inside the co-rotation radius for 
many BHCs at 10^^ erg/s. These would auto- 
matically have hard spectra at this luminos- 
ity 

Three separate regions contribute to the 
spectra of NSs. These are the surface and 
co-rotating magnetosphere, the inner disk- 
magnetopause, and the outer disk. It is ac- 
cepted that blackbody and/or thermal brchms 
spectra of /cTf,f, ~ 2 kcv arises from the re- 
gion near the surface. Hard power- law ra- 
diations seem to come from the inner disk- 
magnetopause region, and soft excess themal 
radiations kTijb ~ 0.1 - 0.6 kev originate from 
the outer disk. Not only are these spectrally 
different regions, they have been spatially re- 
solved in dip observations. Morley et al. 
(1998) have analyzed deep (100%) dips in the 
high state spectrum of XB 1916-053 to show 
that kTfjf, = 1.95 keV blackbody emissions 
arise from a rapidly covered central NS while 
a strong power-law spectrum beyond 100 kev 
with photon index of 1.75 was produced by a 
larger emission region. Observations of hard, 
low states show dominant power-law spectra 
along with a low temperature soft excess that 
produces a few percent of the luminosity. 



Exactly the same spectral features and 
some of the same spatial features have been 
observed for BHCs. For BHCs, however, the 
strong soft radiations are attributed to the 
disk and are usually represented by MBB ra- 
diation functions. It is a stroke of good for- 
tune that MBB functions represent the cen- 
tral soft sources of both NSs and BHCs, but 
certainly no proof of an event horizon. A sim- 
ilar analysis of high state dips of the BHC 
GRO J1655-40 (Kuulkers et al. 1997) re- 
vealed a bright, but unresolved central region 
that produced a power-law and a MBB of 1.1 
kev temperature. The entire x-ray emitting 
region was found to be smaller than ~ 460 
km. The disk diameter at the co-rotation ra- 
dius found here for this source would be 550 
km. The inner disk would have to be smaller 
than this for a NS to produce the high state. 
In a similar way, the times of ingress and 
egress for dips of Cygnus X-1 (Kitamoto et 
al. 1984) constrain the size of the region of 
i^jggin of its low state hard spectrum to be of 
order 4000 km for luminosities for which the 
disk inner diameter would be near 3400 km 
(Misra & Meha 1997). 

Thus the inner disk-magnetopause region 
is a likely site of hard photon production. 
This would be consistent with the large solid 
angles relative to the disk seen for reflected 
x-rays as disk radii expand in decline. Zhang, 
Yu & Zhang (1998) have suggested that flow 
reversal at the magnetopause could produce 
a hard power-law spectrum by the same bulk 
Comptonization mechanism examined by Kluz- 
niak & Wilson (1991), Hanawa (1991), and 
Walker (1992) for flow into a boundary layer. 
Whether or not flow reversal bulk Comp- 
tonization or magnetospherically driven re- 
turn flows (Arons et al. 1984) are related to 
the hard spectrum, the energy requirements 
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for the power-law strongly suggests that it 
would have to originate in the inner disk. 

For both Cygnus X-1 (Balucinska-Church 
et al. 1997) and GRO J1655-40, (Kuulkers 
et al. 1997) a soft excess was found to com- 
prise ~ 6-8% of disk emissions. It was mod- 
eled blackbody in both cases. Tempera- 
tures of 0.13 keV and 0.60 keV, respectively, 
were found for the soft excesses of the two 
sources. Balucinska-Church et al (1995) iden- 
tified the soft excess as disk emissions. The 
relatively stronger presence of the soft excess 
in the deeper dips of GRO J 1655-40 confirms 
that they must arise from the outer disk. It is 
unlikely that all of the gravitational binding 
energy available from the outer disk is used in 
producing a soft excess luminosity . The soft 
excess is quite small compared to the power- 
law which also originates in the disk. However 
if the fraction of binding energy that produces 
the soft excess is roughly independent of ra- 
dius, the disk temperature will still scale as 
T ~ 7--3/4^ ^Yic temperature of 0.13 keV for 
Cygnus X-1 at 1700 km scales up to 0.58 keV 
at the 230 km of GRO J1655-40. This pro- 
vides additional confirmation of the outer disk 
as the source of much of the soft excess. 

An ultrasoft spcctrTim was first recognized 
for BHCs by White k Marshall (1984) for 
LMC X-1 and LMC X-3. Analysis of Rossi 
observations (Wilms et al. 1998) of these per- 
sistent, high state sources have shown that 
their spectra can be fitted with a MBB (or 
blackbody as well for LMC X-1) of temper- 
ature, kTiih = 1-0 and 1.25 keV, respectively 
for these sources, along with a power law. A 
5-9 Mq Yilmaz NS emitting as a blackbody 
would have a surface redshift of z ~ 0.5 — 2.4 
(Fig. 1, Table 2). A local surface tempera- 
ture of 2.3 keV would produce a distantly ob- 
served apparent temperature of T/(l-|-z) = 



2 keV for a low mass NS and 1.3 - 0.7 keV 
for the BHC. Thus bright, ultrasoft peaks are 
well explained for massive NSs. 

We conclude that all of the distinguish- 
ing characteristics of BHCs are encompassed 
by the massive NS hypothesis. If a suitable 
model for NSs were available, we might apply 
it with some confidence to the BHCs. Un- 
fortunately, none is available. Blackbody or 
thermal brehmsstrahlung might be expected 
from near the star surfaces. Such simple mod- 
els (as well as apparently inappropriate high 
temperature MBBs) do a fair job of repre- 
senting the soft spectral components of both 
BHCs and low mass NSs. What is lacking is 
an understanding of the power-law emissions. 
They so dominate the low state spectra, leav- 
ing only small soft excesses, that an optically 
thick, viscosity dominated disk (with its as- 
sociated lower temperature MBB spectrum) 
may not be an adequate starting point for the 
hard x-ray spectrum. It is not certain that an 
optically thick disk at the temperature of the 
soft excess could generate a sufficiently ener- 
gized corona to produce the power-law com- 
ponent. 

The low state disk may be optically thin 
for x-rays and/or advective, i.e., not domi- 
nated by viscosity. An advective accretion 
fiow (ADAF) model was proposed (Narayan 
et al 1997) to explain the hard quiescent ra- 
diations of BHCs. In this model, the disk 
is optically thin and very hot, and accretion 
rates above 10^^^^^ g/s are needed to produce 
quiescent luminosity. According to Eq. 14, 
Atoll sources with magnetic fields of ~ 2a;10^ 
G would have magnetospheric radii of 50 - 
100 km for these accretion rates (or less, ac- 
cepting Dong Lai's suggestion). With spins 
of ~ 350 Hz they would produce ~ 5x10'^'*^'^^ 
erg/s for these accretion rates without any 



22 



surface contributions and much more unless 
propeller effects block surface emissions com- 
pletely. Since quiescent luminosities of Atolls 
are far below these levels, the high accretion 
rate ADAF model is simply inapplicable to 
quiescent NSs. Recall that ADAF models 
with transitions to the conventional optically 
thick, geometrically thin, disk (Esin, McClin- 
tock & Narayan 1997) at radii beyond lO^Rg 
have been shown (Zycki, Done & Smith 1997) 
to fail for BHCs as well. Finally, the presence 
of a hard spectrum for lEl 740. 7-2942 for an 
inner disk radius of only 225 km shows that 
the ADAF model simply cannot be taken se- 
riously. While there seems to be little to rec- 
ommend pure ADAF models, a combination 
of partially advective and magnet ospheric ef- 
fects cannot be ruled out. An advective disk 
with a magnetospherically driven return flow 
(Arons et al. 1984) might have some merit. 
Advective flows can have positive energy for 
the return. 

The nature of the spectrum that would 
arise from within the magnetosphere for lu- 
minosities above Lc is also unclear. The pro- 
duction of soft photons via cyclotron emission 
would be difficult to quantify. The spectrum 
of photons arising from the surface might 
be blackbody or thermal brehmsstrahlung. 
Psaltis, Lamb & Miller (1995) have calculated 
spectra and successfully reproduced the color- 
color diagram for a Z-source with an accretion 
disk-central corona-magnetosphere model. It 
possible that this model might be extended 
to encompass the hard spectrum of the pro- 
peller regime, for which the hot corona might 
be located at the magnetopause. The hard 
spectrum of SAX J1808. 4-3658 illustrates an 
additional difficulty for spectral models. The 
emergent observed spectrum may be strongly 
affected by inchnation effects. 



Black hole models must necessarily use an 
optically thick disk to produce soft peaks and 
some Comptonization of soft photons to pro- 
duce power-laws. Although current models 
with central black holes and accretions disks 
have had some success in representing static 
spectra for some states of the BHCs, most 
models rely on features such hot central 
corona introduced in a purely ad hoc fash- 
ion. It remains to be shown rigorously that 
a hot central corona would develop for Kep- 
lerian flow into an event horizon. The cur- 
rent "sombrero" models (Poutanen 1998) do 
a good job of representing some BHC spectra 
but require that the accretion disk intrusion 
into the central corona be "just so" to let just 
the right number of seed photons be Compton 
scattered to produce the hard spectrum. It is 
not clear that such a model can account for 
the timing signatures within the spectrum. 

6. Active Galactic Nuclei 

It must be acknowledged that active galac- 
tic nuclei pose an interesting challenge for 
any alternative theory of gravity. Spectral 
state switches have been observed for them 
(Mannheim et al 1995). They likely do not 
possess star cluster cores (Maoz 1995). Mas- 
sive single objects with accretion disks might 
be possible but they would be most unusual 
astronomical objects. It can be readily seen 
from Eqs. 11 and 12 that a mean density of 
~ Ig/cm^ could correspond to a very com- 
pact object of ~ IO^Mq that would not nec- 
essarily be a black hole, even in SM. Given 
the extremely high efficiency of conversion of 
gravitational mass-energy to luminosity, there 
should be no problem modeling AGN levels 
of radiation. But the constant proper den- 
sity model is wholly inadequate. Realisti- 
cally, an object of such low density might need 
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to be supported internally by radiation pres- 
sure (but see Graber 1998) and would likely 
have large density gradients and strong mag- 
netic stresses. This is clearly a topic for an- 
other time, but one that should be pursued. 
Black holes would certainly be simpler, but 
we should confirm them among the galac- 
tic BHCs, where we can strongly constrain 
masses and radii, if at all possible. The is- 
sue is clearly that of the existence or non- 
existence of an event horizon. YM objects 
with u{R) > 1/2 are more compact than 
black holes of the same mass and might well 
be of low luminosity with their surfaces inside 
the photon orbit. 

7. Conclusions 

The characteristics of BHCs appear to be 
compatible with a massive NS model. The 
signatures of this model are strong magnetic 
fields, large surface binding energies and large 
redshifts of surface emissions. The field can 
produce power-law and gamma emissions, jets, 
spectral state switches and cyclotron emis- 
sions and the surfaces can produce ultrasoft 
peaks and large rms flickering. The magnetic 
field strengths and spins proposed here for 
BHCs are reasonably near the spin-up trend 
line from radio pulsars to millisecond pulsars. 
Where co-rotation radii can be determined 
from spectral fit parameters, the correspond- 
ing spin rates and subsequently found mag- 
netic field strengths imply quiescent luminosi- 
ties in agreement with observations. Due to 
their slower spins, the magnctosphcric contri- 
butions to the quiescent luminosities of the 
BHCs appear to be somehat smaller than 
those of the low mass Atoll and Z NSs. The 
work of Z^ycki, Done and Smith (1998a) de- 
scribes a spectral state switch in exact ac- 
cord with the magnetic propeller mechanism 



cutting off the fiow to the surface with lit- 
tle change of inner disk radius. Cir X-1, the 
NS most like a BHC, has been found to have 
spin and magnetic field similar to many of the 
BHCs. The ubiquitous slow ~ 6 Hz QPOs 
found in high states of both BHCs and NSs 
(Miyamoto, Kimura & Kitamoto 1991, Mak- 
ishima et al. 1986). can be explained as accre- 
tion disk fiow and opacity oscillations driven 
by surface radiations (Fortner, Lamb & Miller 
1989). It is not known how a black hole might 
produce them. 

A clear implication of the massive NS model 
is that hard spectral features arise from the 
inner disk (magnetopause?) while soft spec- 
tral features can arise from both the disk and 
the surface. The strong refiected x-rays seen 
in hard states with fairly large solid angles 
subtended as the disk radius expands in de- 
cline, suggests that the magnetopause might 
be a primary site of hard photon production. 

Adoption of the massive NS hypothesis 
would eliminate many difficulties for accretion 
disk theory. Spectral state switches would no 
longer require discontinuous disk behavior for 
transitions between high and low states. An 
additional burden for disk theory is removed 
by no longer requiring disk explanations for 
occasional BHC failures to produce ultrasoft 
peaks at high luminosities. The massive NS 
hypothesis explains these as either magnetic 
axis orientation effects or failure of the in- 
ner disk to penetrate the co-rotation radius. 
The structure of the transition region near the 
magnetopause is likely to be much more com- 
plex than previously implied here. It is likely 
to be involved in a complicated, inclination 
dependent way in the production of the pairs 
of kHz QPOs. 

Three of the BHC cases examined here, in- 
cluding Cyg X-1, seem to have spins below 
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10 Hz and very strong magnetic fields. All 
are above the "death line" for pulsars, while 
V0332+53, the original flickering pulsar, is 
not. Cyg X-1 is a wind accreter that appears 
to be delicately balanced between intermedi- 
ate and low states. Its large co-rotation ra- 
dius permits a large energy release and conse- 
quent large radiation pressure increase when 
mass does reach the star surface. The net 
effect is a somewhat self-regulated luminos- 
ity. It is likely near spin secular equilibrium, 
as are the Zs and Atolls. Aside from their 
spins being revealed during bursts, the Zs and 
Atolls are not x-ray or radio pulsars. Mag- 
netic fields of order 10^° G and higher may 
be sufficient to suppress bursts, however, Cir 
X-1 has exhibited bursts. The similarly weak 
magnetic flelds of GRO J 1655-40, A0620-00 
and GRS 1915+105 suggests that they should 
be carefully monitored for bursts; indeed this 
last source may have already produced some 
(Taam, Chen, & Swank 1997). That the 
stronger-field BHCs are not bursters would 
be expected. The lack of pulses requires near- 
coincident magnetic and spin axes or extreme 
magnetic fields. Such a configuration is likely 
necessary for efficient production of jets and 
may enhance gamma emissions. 

KHz QPOs corresponding to the Keplerian 
frequencies of the inner disk imply that the 
accretion disks of both NSs and BHGs can 
reach to near the innermost marginally sta- 
ble orbit. With a fixed binding energy for 
this orbit it seems clear that we must look 
for the differences between NSs and BHGs in- 
side this orbit. Black hole disk models that 
are differentiated from NS disks without ex- 
plicitly incorporating the differences between 
event horizons and star surfaces should have 
little credibility. 

A tremendous amount of detailed model- 



ing will be necessary to extend the massive NS 
hypothesis beyond this introduction. The Yil- 
maz modifications of general relativity have 
provided an acceptable basis for the explo- 
rations of the differences between surfaces and 
event horizons. The simplicity of the Yilmaz 
theory, both mathematically and conceptu- 
ally, warrants its continued use. Its strong 
metric similarities to general relativity below 
2Mq and the lack of an event horizon be- 
yond 2>Mq will provide for stringent compar- 
isons of black hole vs NS models. The re- 
cent discovery of NSs of approximately 2 Mq 
(Zhang et al. 1998) should permit detailed 
comparisons of the two theories for NSs in 
a range where differences might first appear. 
For example, the amplitudes of QPOs gener- 
ated from bursts might be expected to differ 
due to differences of gravitational curvature 
of photon trajectories. Comparisons of this 
sort might be comphcated for masses large 
enough for the surfaces to be inside the pho- 
ton orbit. The important point, however, is 
that the Yilmaz metric provides a sharp tool 
that should eventually decide the question of 
the existence of event horizons. 
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Appendix 1. Particle Mechanics 

The mechanics of a particle orbit can be 
examined very simply with the aid of the 
energy-momentum four- vector. The magni- 
tude of this vector, given by g'-'PiPj, is mo^c^ 
where mg is the rest mass of the particle. For 
a particle in an equatorial orbit (^^ = vr, p^ = 
0) about an object of gravitational mass M in 
the SM, one obtains: 



^2 



[I - 2uy 



- (1 - 2u)pr' - ^ = mo'c' (1) 



Here po = E/c, where E is the particle energy 
and u = GM/c^r is the gravitational poten- 
tial at distance r from the center of mass M. 
P0, the particle angular momentum, is a con- 
stant of the motion. Further, by defining a = 
(cp^/GMmo) and rearranging, the equation 
above becomes: 



2u) 



2 Pr 



;1-2m)(1 + aV) 
(2) 

For suitably small energy, bound orbits occur. 
Turning points for which Pr = can be found 
by examining the effective potential, which 
consists of all terms to the right of E^/mo^c^ 
above. At minima of the effective potential 
we find 

= (3) 
and with p^ = 0, we get 



E = moc 



2 (1- 
V(l - 3^) 



(4) 



The innermost marginally stable orbit is found 
by setting the first two derivatives of the effec- 
tive potential with respect to u to zero. This 
gives the well-known results, u = 1/6 and a? 
= 12. With this value of u, E = ^(8/9) moc^ 



Substituting this value for E, the radial mo- 
tion equation for a subsequent geodesic fall to 
the star surface is: 
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Pr 
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Equating p^ to mv^, using m = E/c^ and rec- 
ognizing (1 - 2u)vr as the proper velocity at 
the position where the potential is u, the ra- 
dial velocity relative to a distant observer is 
seen to be 



(6) 



Analogous treatment of the orbit equation 
in the YM yields: 



exp(— 4m) 



Pr 



E^ 



-exp(— 2m)(1 + a u exp 
(7) 



Setting the first two derivatives of the effec- 
tive potential to zero in this case produces 
coupled equations which have a solution for 
« = (3 - v^)/4 ^ 0.191 with = 12.4. Cir- 
cular orbits occur for 



exp (2m) 
u — 2u^ 



and 



E — moc^ exp(— m) 



l-u 



(8) 



(9) 



. 1-2m 

A particle in geodesic infall starting with p^ = 
at the marginally stable orbit would arrive 
at the star with velocity components: 

v'r = c[l - 1.12exp(-2M)(l + 12.4^2 exp(-2n)]^/2 

(10) 

t;^ = 3.73cm exp(-2M) (11) 

The radial velocity is supersonic for disk tem- 
peratures of as much as 1 keV for u = 0.2; i.e., 
for very httle change of potential relative to 
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that of the marginally stable orbit. For u = 
0.25, corresponding to the surface of a star 
of ~ 3 Mq, the radial speed for a Keplerian 
fall from the marginally stable orbit would be 
about 10^ cm/s. 

The Keplerian frequency of the marginally 
stable orbit can be found from the relations: 

= GMmoa and 



27r dt 

f = ^exp(-2M) to be: 



2nm^)C dt ' 



VGMexp{-2u) 
27rr3/2 J(l - u) 



(12) 



In either metric there is an unstable circu- 
lar photon orbit. In the YM this lies outside 
the star surface for u(R) > 1/2, or for m^ > 
6 Mq. For black holes the unstable photon 
orbit occurs for radius r = 1.5 R^. Although 
radially directed photons can always escape 
in the YM, those with too much orbital an- 
gular momentum can be trapped inside the 
unstable photon orbit. 

2. Compact Objects 

The constant proper density object is one of 
constant intrinsic local density. Since there 
are cases of interest for extensions of this 
model to variable local density, the appropri- 
ate equations will be set up for a more general 
case. Consider a spherically symmetric den- 
sity distribution poh(r<R) with po ^ constant 
and h(r) a distribution function for local den- 
sity. The condition h(0) = 1 is imposed with- 
out loss of generality, and all mass is confined 
within radius R relative to a distant observer. 
In this model, the proper density, p(r)g(r) 
is assumed to be given by poh(r). Thus it 
will be assumed that p(r) = poh{r)/g{r). As 
shown by Clapp (Clapp 1973), u(r) can be 
constructed by adding the potentials of suc- 
cessive shells of matter. Inside each shell the 



contribution to u(r) is constant and outside it 
declines as 1/r. For a general radius r, with 
r < R, contributions to the potential must be 
separated into those arising from shells of ra- 
dius r' inside r and those from shells of radius 
r' outside r. The gravitational mass within 
each shell is Aixr''^ pQh{r')dr' / g{r'). (For the 
free baryon mass substitute g(r'Y^'^ .) The 
central potential u(0) is obtained by dividing 
each shell mass increment by r' and integrat- 
ing out to R. Subtracting u(0) from u(r<R) 
leads to Clapp's integral equation: 

u[r < R)—u[0) = 



dr'h{r'){r' — r' /r)/g{r') 
(13) 



It can be shown that this equation is a solu- 
tion of the field equations of either gravita- 
tional theory by substitution of the appropri- 
ate g(r) and other metric components into the 
field equations and solving the Gtt differen- 
tial equation for u(r). By expressing u(r<R) 
and u(0) in terms of g(r) and g(0) and incor- 
porating g(0) into a scale factor for lengths, 
the integral can be recast in the dimensionless 
form: 

yix) = -2 r f{x')hix')ix' -x'^/x)dx' 
Jo 

(14) 

where x = r/(rog(0)^/^), ro = c/^/AirGpo, f(r) 
= g(0)/g(r), N = 1 in the YM and N = 2 in 
the SM Thus these collapsed objects can be 
described with the aid of the characteristic 
radius, ro, and a characteristic mass, Mq = 
c^Tq/G. In the exponential metric: 



y{x) = ln{f{x)) 



and in the SM: 



y{x) = 1 - l/f{x) 



(15) 



(16) 



After substituting the appropriate left mem- 
ber, the integral equation can be numeri- 
cally solved for f(x) by an iterative relaxation 
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method starting from f(x) = 1 for x < 1 
and an asymptotic value for x > 1. Asymp- 
totic solutions for the constant intrinsic den- 
sity model, for which h(r) = 1 are f(x) = 
for the YM, (Clapp 1973) and found 
here for the SM. Physical quantities such as 
gravitational mass rrig, equivalent free baryon 
mass,m5, distantly observed red shift of sur- 
face radiations, z and radius, R can be ex- 
pressed in terms of certain integrals over f (x) . 
Following Clapp (Clapp 1973), these are: 

FAx)^ r hix')fix')x'dx' (17) 

F2ix)^- r hix')fix')x'^dx' (18) 
X Jo 

F3(x)-- r f{x'f'^x'^h{x')dx' (19) 
X Jo 

In the exponential metric there follows: 

mg = MoxF2{x) exp(-Fi(a;)) (20) 

mft = MoxFsix) (21) 
i? = roxexp(-Fi(x)) (22) 
u{0) = F,{x) (23) 
u{R) = F^ix) (24) 
In the SM the same quantities are given by: 

^F,{x) ^25) 

(26) 
(27) 

(28) 

(29) 



rrik = Mq 
R 



ro- 



uiO) 
u{R) 



°(l + 2Fi(x))2 

xF-jjx) 
(l + 2Fi(x))3/2 

x 

l + 2Fi(^)) 

Fijx) 
l + 2Fi{x) 

l + 2Fi(x) 



Numerical solutions of Eq. 14 and the ap- 
propriate integrals for the case h(r) = 1 were 
used to generate the data of Table 2 and Ta- 
ble 3. Plots of mass, radius and redshift for 
surface emissions are shown in Figures 1 and 
2. Choosing h(r) = 1 forces the object to 
be of constant local intrinsic density. To an 
external observer, there would be substantial 
variations of density, but locally everything 
would seem entirely normal. Such is life in 
radically curved space-time. Simple exten- 
sions of this approach to an Oppenheimer - 
VolkofT method for the determination of po 
and h(r) are possible, as well as exponential 
atmosphere models of active galactic nuclei 
(h(r) not equal to 1), but the essential differ- 
ences between black holes, neutron stars and 
nuclear-density compact objects in the expo- 
nential metric are clearly revealed by this sim- 
ple model. 
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Table 1 
Neutron Star Properties 
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^Values calculated in this work. 
^QPO frequency. 
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Table 2 

Quantities in the Yilmaz Metric 
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Note.- 


— For X, f(x^ 


), sec appendix. R 


is in units of 


ro = c 




and mg 


in units of 



Mo = c^tq/G. g is proper surface free-fall acceleration. 
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Table 3 

Quantities in the Schwarzschild Metric 
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0.007 


0.266 


0.027 


0.026 


0.038 


0.099 


0.498 


0.925 


0.028 


0.026 


0.401 


0.071 


0.062 


0.097 


0.172 


0.653 


0.880 


0.050 


0.044 


0.466 


0.111 


0.095 


0.143 


0.226 


0.803 


0.832 


0.073 


0.061 


0.505 


0.149 


0.122 


0.185 


0.277 


0.918 


0.793 


0.090 


0.073 


0.524 


0.178 


0.140 


0.214 


0.314 


^.923 


0.792 


0.091 


0.074 


0.525 


0.180 


0.141 


0.215 


0.316 


1.003 


0.765 


0.102 


0.081 


0.534 


0.199 


0.152 


0.234 


0.341 


1.153 


0.716 


0.122 


0.093 


0.545 


0.231 


0.170 


0.264 


0.385 


1.353 


0.655 


0.143 


0.104 


0.550 


0.269 


0.190 


0.297 


0.438 


1.513 


0.609 


0.157 


0.111 


0.550 


0.295 


0.202 


0.318 


0.475 


1.653 


0.573 


0.168 


0.115 


0.549 


0.314 


0.210 


0.334 


0.503 


1.853 


0.526 


0.180 


0.120 


0.546 


0.336 


0.220 


0.353 


0.539 


2.053 


0.484 


0.190 


0.123 


0.542 


0.354 


0.227 


0.368 


0.569 


2.353 


0.430 


0.201 


0.126 


0.536 


0.375 


0.236 


0.386 


0.604 


2.753 


0.373 


noil 


n 1 0S 


n i^oQ 
u.ozy 


u.oyo 


n 0/10 


n /in/i 


u.uoy 


3.053 


0.338 


0.216 


0.129 


0.524 


0.403 


0.246 


0.414 


0.658 


3.553 


0.291 


0.222 


0.129 


0.518 


0.413 


0.250 


0.427 


0.680 


4.053 


0.255 


0.226 


0.129 


0.514 


0.419 


0.252 


0.437 


0.695 


4.998 


0.206 


0.231 


0.129 


0.508 


0.424 


0.253 


0.449 


0.710 


6.050 


0.168 


0.233 


0.129 


0.504 


0.429 


0.255 


0.458 


0.725 


8.050 


0.125 


0.235 


0.128 


0.501 


0.427 


0.254 


0.469 


0.725 


10.05 


0.100 


0.235 


0.127 


0.500 


0.425 


0.254 


0.475 


0.722 


15.05 


0.066 


0.236 


0.126 


0.500 


0.421 


0.252 


0.483 


0.717 


20.00 


0.050 


0.236 


0.126 


0.500 


0.419 


0.252 


0.488 


0.713 


^Entries for 1.4 Mq 


with po = 


1.35a;10" 


g/cw?. 











Note. — For x, f(x), see appendix. R is in units of rg = c/\/47rGpo, and in units of 
Mq = c^Tq/G. g is proper surface free-fall acceleration. 
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